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摘 要 : 在 西 营 河 九条 岭 站 以 上 流域 内 ,利用 2011 一 2018 年 格 点 降水 和 气温 数据 集 的 逐日 资料 ,通过 降水 径流 成 


因 关系 探索 出 格 点 降水 异常 点 修正 方法 ,并 将 融雪 模块 与 三 水 源 新 安江 模型 耦合 ,研究 流域 日 径流 模拟 过 程 。 采 


] 确 定性 系数 和 相对 误差 两 个 评估 准则 进行 模拟 效果 的 记 


F 价 ,将 降水 修正 后 以 及 融雪 耦合 后 的 模拟 效果 与 未 改进 


的 模型 模拟 效果 进行 比较 ,分 析 融 雪 模 块 耦合 后 的 新 安江 模型 在 西 营 河 流域 的 适用 性 。 结 果 表 明 :参照 降水 径流 
成 因 关 系 进行 的 格 点 降水 修正 提高 了 径流 模拟 的 精度 ,75% 的 年 份 模拟 径流 的 确定 性 系数 有 所 提高 。 耦 合 融雪 的 


新 安江 模型 在 研究 区 模拟 较 好 ,75% 以 上 年 份 的 模拟 比 未 改进 时 模拟 效果 好 ,87% 以 上 年 份 的 模拟 确定 性 系数 达 


0.6 以 上 上。 耦合 融雪 的 新 安江 模型 为 西 营 河 流域 融雪 径流 的 预报 预警 提供 决策 参考 。 
关键 词 : 西 营 河 流域 ; 格 点 降水 ; 融雪 模块 ;改进 的 新 安江 模型 


中 国 西北 地 区 的 河流 多 发 源 于 山区 ,冰雪 融 水 
是 河流 春季 的 重要 补给 来 源 沾 。 西 营 河 位 于 甘肃 省 
境内 ,其 径流 主要 来 源 于 祁连山 东 段 融雪 和 山区 降 
水 。 在 西 营 河流 域 ,多 数 研究 侧重 于 西 营 灌 区 的 生 
态 及 水 利 建设 ,进行 水 文 模拟 极 少 ,迫切 需要 探索 
适用 于 该 区 域 的 水 文 模型 。 在 该 流域 范围 内 使 用 
水 文 模型 进行 模拟 时 ,需要 考虑 春季 融雪 的 影响 。 
Wu 等 通过 研究 长 序列 融雪 径流 过 程 ,完成 了 融雪 
组 分 的 评估 ,指出 春 夏季 融雪 径流 比例 随 气 候 变 暖 
而 有 所 改变 ,融雪 过 程 会 影响 高 寒山 区 水 文 过 程 。 
李 志 龙 2 在 研究 资料 缺乏 的 寒 区 流域 水 文 模 拟 中 ， 
将 改进 的 新 安江 模型 作为 混合 模型 的 一 部 分 ,进行 
了 融雪 径流 的 模拟 。 目 前 ,关于 融雪 径流 预测 已 经 
有 应 用 较 广 的 模型 ,如 SRM 模型 六 ;此 外 还 有 包含 融 
雪 模 块 的 模型 ,如 SWAT 模 型 5 MIKE SHE 模型 “5 
等 。 然 而 ,SRM 模型 主要 对 于 春季 融雪 模拟 较 好 ， 
多 数 只 能 应 用 于 年 内 短期 的 径流 预报 ”。SWAT 模 
型 与 SHE 模型 又 较为 复杂 ,输入 的 数据 需要 包含 大 
量 信息 (包括 数字 高 程 数 据 .土地 利用 数据 、 积 雪 
IERIE .土壤 特征 .降水 气象 数据 等 ) ,在 西 营 河 这 


少 ,代表 性 不 强 。 故 而 需要 探索 更 加 适用 于 西 营 河 
流域 的 融雪 径流 模型 。 

新 安江 模型 * 作 为 国内 应 用 较 广 的 模型 之 一 ， 
其 模型 结构 清晰 、 层 次 分 明 ,在 实际 应 用 中 得 到 不 
断 地 完善 ,可 以 通过 改进 的 模型 进行 数据 缺乏 地 区 
的 水 文 模拟 。 刘 金涛 等 "总 结 了 新 安江 模型 发 展 
的 特点 ,主要 是 架构 形式 趋 于 多 样 性 、 产 汇流 模型 
不 断 物理 化 .应 用 范围 不 断 延 伸 。Lyu Se a 
子 群 算法 和 集合 卡尔 曼 滤波 数据 同化 的 3 种 变化 ， 
对 罗江 和 资料 匮乏 地 区 进行 模型 参数 批量 佑 计 。 
韩元 元 等 "将 变动 态 存储 系数 法 结合 新 安江 模型 
进行 河道 洪水 演算 。 以 上 研究 说 明 在 模型 驱动 数 
据 及 方法 上 均 可 与 新 安江 模型 进行 配合 ,从 而 提高 
模型 的 应 用 范围 及 模拟 效果 。Alazzy "采用 不 确定 
性 佑 计 (GLUE ) 方 法 ,对 中 国 新 安江 模型 估算 月 径 
流 的 不 确定 性 进行 检验 和 预测 ,发 现 采 用 Nash 效 率 
系数 评估 的 新 安江 模型 的 参数 不 确定 性 更 小 。 据 
此 ,本 文 率 定 参数 时 ,主要 以 Nash 效率 系数 作为 率 
定 参 数 时 的 目标 函数 。 邓 元 倩 等 ”、 俘 成 国 等 "分 
别 将 新 安江 模型 应 用 于 尘 河 流域. 新 疆 乌鲁木齐 


种 数据 较为 缺乏 的 流域 ,数据 搜集 困难 且 数 量 较 
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扩大 到 西北 地 区 ,不 青 仅 仅 局 限于 湿润 地 区 。 赵 兰 
兰 等 “根据 干旱 地 区 实际 情况 ,将 超 渗 产 流 模式 加 
入 新 安江 模型 在 拒 马 河流 域 进行 了 暴雨 洪水 预报 
及 对 比分 析 , 在 新 安江 模型 各 层次 内 部 使 用 的 方法 
可 根据 流域 的 地 理 位 置 等 条 件 进 行 适 当 改 进 。 至 
此 , 超 渗 产 流 与 新 安江 模型 的 融合 ,使 新 安江 模型 
突破 了 原 有 的 固化 模式 ,从 最 根本 的 产 流 层次 进行 
改进 ,加 强 了 新 安江 模型 在 全 国 多 地 的 适用 性 。 美 
开 芳 "最 早 将 融雪 径流 模拟 加 入 新 安江 模型 应 用 
于 切 德 克 流域 ,随后 田 龙 "等 人 在 该 版 本 基础 上 考 
虑 能 量 概念 进行 高 寒山 区 融雪 径流 模拟 。 然 而 ,该 
模型 仍 需 对 太阳 辐射 .地形 影响 、 积 雪 面 积 等 进行 
量 分 析 ,在 资料 缺乏 的 山区 应 用 较 难 。 本 文选 取 
应 用 广泛 .层次 结构 化 清晰 的 新 安江 模型 ,在 其 中 
加 入 仅 以 气温 为 主要 阔 值 的 融雪 模块 ,从 而 扩大 新 
安江 模型 的 适用 范围 。 

在 资料 匮乏 地 区 ,如 西 营 河 九条 岭 站 以 上 流域 
无 降水 观测 , 且 流 域内 地 形 差异 大 ,九条 岭 水 文 站 
降水 无 法 代表 全 流域 情况 。 较 为 准确 的 降水 驱动 
数据 对 水 文 模拟 具有 重要 意义 ,直接 影响 水 文 模拟 
精度 。 在 模拟 时 ,需要 选取 最 适合 的 驱动 数据 。 孙 
铭 悦 ”等 在 进行 水 文 模型 模拟 时 ,比较 了 站 点 数据 
和 格 点 降水 数据 的 模拟 效果 ,得 出 格 点 降水 比 站 点 
平均 的 降水 数据 更 好 。 本 文 首先 进行 格 点 降水 数 
据 修正 ,随后 将 融雪 模块 与 新 安江 模型 耦合 ,利用 
修正 后 的 降水 数据 驱动 耦合 模型 进行 径流 模拟 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

西 营 河 位 于 武威 市 西南 部 ,是 石 羊 河流 域 的 一 
条 支流 ,上 游 由 宁 昌 河和 水 关 河 组 成 ,发 源 于 祁 连 
山东 段 冷 龙 岭 。 西 营 河 属于 高 寒 半 干旱 半 湿 润 气 
候 区 。 径 流 主要 由 高 山 融雪 以 及 祁连山 区 降水 组 
成 。 西 营 河 多 年 平均 出 山 径流 达 3.184x10* me ,位 列 
石 羊 河 上 游 八 条 河流 之 首 , 占 比 22.3%。 出 山口 附 
近 设 有 九条 岭 水 文 站 ,也 是 该 流域 出 山口 以 上 唯一 
的 水 文 测 站 。 九 条 岭 站 以 上 区 域 植被 良好 ,水 土 流 
失 轻 微 ,此 处 与 上 游 由 于 地 势 的 差异 ,降水 差别 较 
大 "站 。 上 游 年 降水 量 可 达 500~800 mm ,九条 岭 站 年 
降水 量 仅 为 200~400 mm。 流 域 4 一 5 月 多 为 春季 融 
雪 径 流 ,6 一 9 月 径流 为 山区 降水 产生 。 本 文 主要 研 


究 区 为 西 营 河 九条 岭 站 以 上 流域 (图 1) , 集 水 面积 
为 1077 km 。 该 站 是 石 羊 河流 域 的 八大 河流 中 不 受 
人 工 设施 影响 的 水 文 测 站 之 一 ,积累 了 长 序列 水 文 
实测 资料 。 
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图 1 研究 区 域 示 意图 
Fig.1 Schematic diagram of the study area 


1.2 数据 来 源 

本 文 改 进 的 新 安江 模型 输入 的 降水 和 气温 数 
据 采 用 国家 气象 信息 中 心 发 布 的 数据 集 。 降 水 采 
用 中 国 地 面 降水 日 值 0.5°x0.5° 格 点 数据 集 (V 2.0), 
代码 为 SURF_CLL CHN_PRE_DAY_CRID_0.5。 气 
温 采 用 中 国 地 面 气温 日 值 0.5°x0.5° 格 点 数据 集 
(V2.0) ,代码 为 SURF_CLI_CHN_TEM_DAY_GRID_ 
0.5。 此 外 ,九条 岭 水 文 站 2011—2018 E WZ H 3$ 
发 和 逐日 流量 资料 来 源 于 中 华人 民 共 和 国 水 文 年 
鉴 第 10 卷 第 5 册 匡 肃 河 西 地 区 内 陆 河水 文 资 料 。 
1.3 研究 方法 
1.3.1 格 点 降水 数据 修正 方法 ” 格 点 数据 集 主要 由 
中 国 地 面 高 密度 台 站 的 降水 结合 薄 盘 样 条 法 制 
作 , 依 托 于 站 点 数据 ,插值 时 存在 一 定 误差 。 数 据 
集 (V 2.0) 制 作 人 赵 煜 飞 纪 对 该 数据 集 进行 精度 评 
佑 时 指出 ,夏季 误差 较 其 他 季节 较为 明显 。 将 格 点 
降水 数据 用 于 降水 (融雪 ) 径 流 模 拟 时 发 现在 夏季 
会 出 现 少量 径流 异常 值 ,需要 进行 降水 数据 修正 。 
经 试验 ,未 经 修正 的 格 点 降水 ,在 研究 区 每 年 都 会 
出 现 少数 的 异常 情况 ,径流 未 出 现 峰值 ,然而 作为 
径流 主要 成 因 的 降水 却 出 现 较 大 值 , 旦 与 序列 前 后 
日 突变 较 大 ,与 临近 站 点 降水 及 相关 记录 比较 发 现 
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此 日 并 非 又 十 日 且 前 后 日 并 无 较 大 降雨 ,从 而 考虑 
应 是 数据 集 制作 产生 的 误差 导致 。 常 用 的 修正 方 
式 是 与 附近 站 点 序列 进行 比较 。 然 而 ,在 西 营 河流 
域 ,附近 的 气象 站 点 距离 流域 较 远 (图 2), 且 分 布 多 
偏 东北 、 西 北方 向 ,作为 降水 产 流 主要 来 源 的 西南 
部 山区 没有 代表 性 测 站 ,无 法 全 面 展 现 整 个 流域 降 
水 .气温 等 情况 。 
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图 2 西 营 河 流域 邻近 气象 站 点 分 布 
Fig. 2 Distribution of meteorological stations near Xiying 


River Basin 


本 文 探索 性 尝试 通过 年 最 大 降水 和 径流 系数 
进行 约束 ,筛选 出 每 年 的 异常 格 点 降水 值 进行 修 
正 , 从 而 提高 径流 模拟 的 精度 。 经 过 多 年 格 点 降水 
径流 模拟 试验 ,归纳 得 出 以 下 降水 量 和 径流 系数 两 
个 约束 能 检测 出 较 明显 的 格 点 降水 异常 数据 。 同 
时 违背 降水 量 约束 (公式 1a) 和 径流 系数 约束 (公式 


1b) 两 个 约束 条 件 判 断 为 异常 点 : 
0< P<O7P i max (la) 
a>0.3, as Ë (1b) 


式 中 : 民 表 示 出 口 九条 岭 水 文 站 日 尺度 实测 径流 深 
(mm);P 为 逐日 格 点 降水 值 (mm ) ;a 为 径流 系数 , 径 
流 系数 0.3 的 约束 条 件 主要 是 采用 3c 原 则 ,通过 正 
态 分 布置 信 区 间 算 出 , 即 平均 年 径流 系数 为 0.503， 
标准 差 c 为 0.067, 在 小 于 w-3c 即 0.3 时 ,概率 仅 为 
0.27% ,认为 该 测量 值 为 坏 值 。 

尽管 下 垫 面 条 件 和 气候 条 件 的 复杂 性 会 对 不 
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同 自然 地 理 条 件 的 降水 -径流 关系 产生 一 定 影响 ， 
降雨 和 径流 仍 具 有 较 强 的 成 因 关 系 ,可 以 作为 两 个 
序列 的 初步 评定 参考 。 本 研究 识别 出 异常 点 后 , 利 
用 降水 和 径流 相关 关系 进行 修正 ,将 识别 出 的 异常 
点 修正 到 每 年 绘制 的 剔除 异常 点 的 降水 -径流 相关 
关系 图 上 ,具体 流程 如 图 3 所 示 。 

某 些 年 份 单 日 格 点 降水 最 大 值 达到 40 mm 甚至 
更 大 时 , 阁 此 点 为 异常 点 ,会 影响 降水 量 约束 的 异 
第 点 挑选 ,检测 时 遇 此 情况 ,需要 以 年 第 二 大 降水 
值 作为 降水 量 约束 的 基础 值 ,再 乘 以 试验 率 定 出 的 
比例 系数 0.7。 
1.3.2 融雪 模块 与 新 安江 模型 耦合 “本文 在 新 安江 
模型 基础 上 ,增加 了 融雪 模块 ,主要 包括 融雪 模块 、 
蒸 散 发 计算 、 匡 满 产 流 计算 、 三 水 源 划 分 和 汇流 
计算 


(1) 融雪 模块 

改进 后 的 新 安江 模型 在 第 一 层次 前 ,加 入 融雪 
模块 ,主要 采用 度 日 因子 法 进行 融雪 计算 。 度 日 法 
假设 融雪 量 与 气温 为 显著 线性 关系 ,其 计算 公式 
如 下 : 


_{dx(T-T,), T>T, 
0 fet, 
AEP: MW fit 3 7k Emm. d’); d HEH NF 
(mm C- d); To Ay fa a FR BE CC); THA 
(Co 
在 融雪 模块 中 ,主要 输入 数据 为 气温 和 降水 。 
积 雪 消 融 的 能 量 主 要 来 源 于 净 辐 射 和 感 热 通 量 ,其 
中 , 净 辐 射 包含 长 波 辐射 和 短波 辐射 。 冰 雪 消 融 以 
近 地 表 层 的 长 波 辐射 和 感 热 通 量 原因 为 主 ,气温 叉 
是 地 表 长 波 辐射 的 主要 影响 因素 ,从 而 本 方法 将 能 
量 的 影响 直接 由 较 易 获取 的 气温 作为 模型 的 驱 
动 。 加 入 融雪 模块 后 ,将 融雪 水 量 和 降雨 量 综 合 的 
液态 水 作为 后 续 模 型 的 输入 。 
(2) 蒸 散发 计算 
模型 采用 三 层 蒸发 模式 ,其 中 , 葵 散 发 能 力 使 
用 燕 发 折算 系数 与 九条 岭 水 文 站 观测 值 的 乘积 。 
EP= KCx Epan (3) 
UH EP AAA ET (mm); KC 为 蒸 散发 折算 系 
数 ;Epan 为 蒸发 阵 观 测 蒸发 mm)。 运 用 的 主要 思 
想 是 土壤 蒸发 能 力 与 相同 气象 条 件 下 的 水 面 蒸发 
一 般 呈 线性 关系 。 
(3) AW ht 


M (2) 
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径流 呈现 负 相 关 


年 最 大 降水 值 与 对 应 实测 


将 年 降水 第 二 大 值 作为 Prnso 后 续 将 此 
日 最 大 降水 值 修正 到 相关 图 


用 降水 与 实测 径流 深 绘制 
降水 -径流 的 相关 关系 图 


绘制 y=0.7Piewox 和 y=x/0.3 线 


将 偏离 曲线 较 大 的 异常 点 选 出 ， 
区 域 : {x,y)ly>0.7Pormers Y>X/0.3} 


用 不 包含 异常 点 的 格 点 降水 与 实测 径流 
重新 绘制 降水 -径流 的 相关 关系 图 


| 


| 将 找 出 的 异常 点 根据 实测 径流 深 修正 到 


新 绘制 的 降水 -径流 相关 关系 图 


是 否 仍 有 偏离 
曲线 较 大 点 


修正 后 的 格 点 降水 数据 作为 
水 文 模型 输入 


图 3 格 点 降水 异常 点 识别 与 修正 流程 


Fig. 3 Grid precipitation anomaly identification and correction process 


aT aR ed ae BF 5 7 E 
流 模式 。 概 化 为 达到 田间 持 水 量 后 , 才 开始 产 流 。 
由 于 流域 各 处 蓄 水 能 力 不 同 ,降水 与 初始 包 气 带 蓄 
水 量 之 和 小 于 流域 包 气 带 蕾 水 容量 最 大 值 时 ,为 流 
域 部 分 面积 产 流 。 同 理 , 降 水 与 初始 包 气 带 蓄 水 量 
之 和 大 于 流域 包 气 带 蓄 水 容量 最 大 值 时 为 全 流域 
产 流 。 具 体 计算 公式 如 下 : 


R=PE+W -wm+wm(1- PEHA) A+PE<WMM 

0 WMM ) ° 
R=PE+W,-WM ~ A+PE>WMM 
(4) 


SUE AA ii SA a Le BK PTT AS os 
值 (mm) ; WMM Yin BEES eK EKRE (mm) ; WM 
Frit EY La AKRE (mm ) ; BC WTS 
包 气 带 蓄 水 容量 分 布 的 不 均匀 性 ,8 越 小 表示 越 均 
匀 ;R 为 产 流量 (mm) ;WW 为 流域 蓄 水 量 (mm)。 值 得 
注意 的 是 ,PE 中 的 P 指 的 是 融雪 计算 后 的 液态 水 量 
(mm). 

(4) 三 水 源 划分 

产 流量 的 划分 延续 传统 新 安江 模型 的 思路 ,分 
为 地 表 径 流 . 培 中 流 和 地 下 径流 。 考 虑 产 流 面 积 的 


影响 和 下 扑面 条 件 , 采 用 与 上 一 模块 包 气 带 匡 水 容 
量 分 布 不 均 类 似 的 处 理 , 仍 将 产 流 面积 上 的 自由 水 
容量 分 布 用 抛物 线 近似 。 其 中 ,自由 水 是 毛管 水 和 
重力 水 的 合 称 ,是 在 土壤 中 受 力 的 作用 可 以 自由 流 
动 的 水 。 
(5) 汇流 计算 
地 面 径流 直接 流入 河 网 ,计算 公式 如 下 : 
QS) = RS@ x U (5) 
式 中 ;05 为 地 面 径流 汇流 (ms ;RS 为 地 面 径 流 
(mm) ;UU 为 单位 转换 系数 ,为 常数 ;i 指 日 数 ,时 间 尺 
度 为 日 尺度 。 
壤 中 流 汇 流 采 用 线性 水 库 法 ,计算 公式 如 下 : 
QI(i)= CIx QIG-1)+(1- CI) XRIGXU (6) 
P OLN REPEL Em’ +s"), CLARE PLA IB A 
数 ,RI 为 壤 中 流 径 流 (mm) ,UU 为 单位 转换 系数 ,为 
常数 。 
地 下 径流 汇流 仍 采 用 线性 水 库 法 ,采用 公式 
如 下 : 
QG = CGX OGG- 1)+(1-CG)xXRGWXU (7) 
式 中 :0C 为 地 下 径流 汇流 (ms ;CC 为 地 下 水 消 
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退 系数 ;RG 为 地 下 水 径流 (mm) ;U0 为 单位 转换 系 
数 ,为 常数 。 

1.3.3 评估 准则 根据 水 文 预 报 评价 标准 ,采用 相 
对 误差 RE 和 确定 性 系数 DC 进行 评价 。 二 者 公式 
如 下 : 


sm | x 100% (8) 


2 Oi Un) 
yO. ~ Oi 
式 中 : Qu 为 九条 岭 站 实测 流量 (ms ); 0 为 模 
型 模拟 的 流量 值 (ms ') ;RE 为 相对 误差 ;DC 为 确 
定性 系数 。 


2 结果 与 分 析 


2.1 格 点 降水 数据 修正 结果 分 析 

根据 流域 在 降水 网 格 的 位 置 ,提取 流域 内 格 点 
降水 。 经 过 格 点 降水 异常 点 识别 与 修正 方法 进行 
修正 (图 3)。 在 约束 条 件 所 构成 的 区 域内 寻找 异常 
点 ,每 年 需要 修正 的 格 点 并 不 多 ,但 这 些 格 点 属于 
偏离 降水 径流 成 因 关 系 曲线 较 大 的 降水 数据 点 。 
由 2018 年 的 格 点 降水 的 修正 可 以 看 出 (图 4),2018 
年 仅 有 两 日 的 格 点 降水 需要 修正 , 即 红色 标注 部 分 
(图 4a) ,将 异常 点 处 的 降水 数据 修正 到 剔除 异常 点 
后 的 降水 -径流 关系 曲线 上 (图 4b) ,这 两 个 异常 点 
的 降水 值 分 别 由 24.3 mm 、21.8 mm 修正 为 0.6 mm 和 
4.3 mm。 修 正 后 , 率 定 期 的 模拟 效果 确定 性 系数 从 
0.59 提升 到 0.63 ,提升 了 0.04。 验 证 期 的 模拟 效果 
确定 性 系数 从 0.69 提升 到 0.72 ,提升 了 0.03( 表 1)。 


DC=1- (9) 


30 (a) 2018 年 格 点 降水 异常 点 识别 
25 L y=R/0.3 


J=0.7Premax=17 


y= 0.181 12+ 
2.3618x-0.2855 
ba 


1 2 3 4 5 6 7 8 
径流 深 /mm 
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2013 年 的 格 点 降水 修正 对 确定 性 系数 的 提升 较为 
明显 ,从 0.23 提升 到 了 0.69。 格 点 降水 修正 对 于 
75% 的 年 份 均 表 现 为 模拟 效果 提升 。 
2.2 融雪 径流 模拟 分 析 

分 别 用 融雪 模型 融合 前 后 及 降水 修正 后 的 新 
安江 模型 进行 模拟 ,模拟 效果 见 表 1。 由 表 1 可 知 ， 
融雪 模块 耦合 后 的 模拟 结果 通常 比 传统 的 新 安江 
模型 模拟 得 更 好 。 可 以 注意 到 ,改进 过 程 中 ,确定 
性 系数 的 提高 和 相对 误差 的 降低 并 不 同步 ,这 是 由 
于 在 参数 优化 时 ,选择 了 确定 性 系数 作为 目标 函 
数 , 其 更 能 反应 模拟 过 程 ,体现 流量 模拟 值 的 变化 
在 总 体 变化 中 的 比例 。 而 相对 误差 更 多 地 评价 模 
拟 总 流量 与 实测 值 之 间 的 关系 ,重点 在 总 量 的 控制 
上 。 除 2013 年 由 于 以 确定 性 系数 为 目标 函数 的 参 
数 优化 使 得 确定 性 系数 大 幅 提 升 ,相对 误差 为 
23% ,超过 了 20% 外 ,其 余年 份 的 各 种 情况 模拟 , 相 
对 误差 均 在 20% 以 下 ,在 后 续 研 究 中 车 能 继续 进行 
多 目标 参数 优化 ,有 望 将 2013 年 的 相对 误差 进一步 
降低 。 

选择 2011 一 2016 年 作为 模型 率 定 期 ,2017 一 
2018 年 作为 模型 验证 期 。 率 定期 参数 优化 后 ,确定 
性 系数 达到 0.63, 此 套 参 数 用 于 验证 期 ,确定 性 系数 
为 0.72。 说 明 此 模型 在 西 营 河流 域 具 有 一 定 适 用 
性 。 经 过 融雪 模块 融合 、 格 点 降水 修正 ,87% 以 上 的 
模拟 年 确定 性 系数 在 0.6 以 上 .62% 以 上 的 年 份 模拟 
达 0.69 以 上 ,模拟 结果 较 好 。 新 安江 模型 平均 确定 
性 系数 为 0.62 ,融雪 耦合 和 格 点 降水 修正 后 平均 确 
定性 系数 为 0.67 ,提高 0.05。 将 变动 态 存 储 系数 法 
用 于 新 安江 模型 河道 洪水 演算 的 研究 在 印 江 河流 


(b) 2018 年 格 点 降水 异常 点 修正 
y=R/0.3 


: 20 G=0.1Piem=17 
15 y= —0.165x+ 
a 2.2396x-0.3174 

全 


1 2 3 4 5 6 7 8 
径流 深 /mm 


注 :红色 点 表示 降水 异常 点 。 


图 4 2018 年 格 点 降水 的 修正 


Fig.4 Grid precipitation correction in 2018 
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表 1 不 同情 况 下 新 安江 模型 模拟 结果 


Tab.1 Simulation results of Xin anjiang model under different conditions 


相对 误差 /9% 确定 性 系数 
实测 径流 深 格 点 降水 修正 格 点 降水 修正 
年 份 降水 量 /mm 融雪 模块 融合 “下 、 ”融雪 模块 融合 CARA 
/mm 新 安江 模型 “融雪 使 蕊 融合 后 的 改进 新 安 “。， 新 安江 异型 ”融雪 模块 融合 “后 的 改进 新 安 
后 新 安江 模型 后 新 安江 模型 
oe 江 模型 模型 
2011 639.7 289.61 14.0 7.8 16.1 0.68 0.64 0.69 
2012 573.2 330.76 19.2 17.0 19.0 0.69 0.71 0.72 
2013 495.3 233.34 23.0 0.2 23.2 0.69 0.23 0.69 
2014 733.5 284.99 14.7 12.4 12.4 0.43 0.51 0.51 
2015 556.5 273.28 14.3 6.0 15.2 0.56 0.65 0.61 
2016 593.5 324.96 15.6 11.0 19.3 0.58 0.61 0.64 
2017 664.5 341.85 14.7 2.7 17.1 0.62 0.65 0.72 
2018 543.0 317.50 LAL 13.2 19.1 0.70 0.70 0.71 
率 定 期 3591.7 1739.94 3.6 16.2 1.7 0.58 0.59 0.63 
验证 期 1207.5 659.35 4.0 21.0 16.2 0.67 0.69 0.72 


WR FE iy T 0.0182" ;改进 新 安江 模型 在 乌江 独 木 
河流 域 模拟 时 确定 性 系数 仅 提高 0.022。 与 类 似 研 
究 相 比 ,本 文 运 用 改进 的 新 安江 模型 在 西 营 河流 域 
的 平均 确定 性 系数 提高 0.05 ,效果 显著 。 

KH 2011—2016 年 进行 参数 优化 模型 率 定 Fl] 
用 2017 一 2018 年 进行 模型 的 验证 ,参数 率 定时 采用 
遗传 算法 优化 参数 ,其 中 以 确定 性 系数 作为 算法 的 
目标 函数 进行 率 定 , 表 2 显示 了 优化 后 的 模型 参数 。 

耦合 融雪 的 新 安江 模型 率 定 期 和 验证 期 模拟 
结果 如 图 5 所 示 , 除 第 二 年 洪峰 模拟 偏 低 外 ,其 余年 


份 洪峰 模拟 效果 均 很 好 。 率 定期 的 多 年 模拟 结 
中 ,确定 性 系数 达到 0.63 , 比 未 进行 格 点 降水 修正 的 
0.59 提高 了 0.04。 在 验证 期 ,可 以 看 到 模型 已 模拟 
出 春季 融雪 径流 ,夏季 的 洪水 部 分 前 期 模拟 较 好 ， 
后 期 模拟 效果 较 低 ,可 能 与 当年 的 格 点 降水 数据 仍 
存在 误 益 有 关 。 验 证 期 耦合 融雪 的 新 安江 模型 经 
降水 修正 后 比 新 安江 模型 模拟 的 确定 性 系数 提升 
0.05。 

在 融雪 径流 的 模拟 过 程 中 , 率 定期 内 ,2011 年 、 
2013 年 .2014 年 .2015 年 的 春季 径流 模拟 很 好 。 在 


表 2 含 融雪 的 新 安江 模型 参数 及 物理 意义 


Tab. 2 The parameters of snowmelt-Xin’ anjiang model and its physical significance 


层次 参数 参数 意义 敏感 性 取 值 范围 参数 取 值 
融雪 模块 To 基础 温度 /%C 敏感 -1.0~2.7 1.13 
d 度 日 因子 /mm.%.d) 不 敏感 0.5~7 0.69 
蒸 散发 计算 KC 流域 蒸 散 发 折算 系数 敏感 <1 0.08 
UM 上 层 张 力 水 容量 /mm 不 敏感 10~50 42.3 
LM 下 层 张力 水 容量 /mm 不 敏感 60~90 63 
C 深层 蒸 散发 折算 系数 不 敏感 0.1~0.2 0.1 
产 流 计算 WMM 流域 平均 张力 水 容量 不 敏感 70~200 117 
B TKI A SAK AS MATT IR 不 敏感 0.1~0.4 0.28 
IM 不 透水 面积 占 全 流域 面积 的 比例 不 敏感 0.01~0.02 0.018 
水 源 划分 SM 表层 自由 水 苹 水 容量 /mm 敏感 10~80 69.4 
EX KJE Arak eek At MATT K 不 敏感 1~1.5 1.01 
KG 表层 自由 水 蕾 水 库 对 地 下 水 的 日 出 流 系数 敏感 <0.7 0.3 
KI 表层 自由 水 琴 水 库 对 壤 中 流 的 日 出 流 系数 不 敏感 0.7-KG 0.4 
汇流 计算 CI 壤 中 流 消退 系数 敏感 0~0.9 0.9 
CG 地 下 水 消退 系数 不 敏感 0.98~0.998 0.998 
CS 河 网 蓄 水 消退 系数 敏感 0.001~1 0.69 
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图 5 模型 率 定 期 和 验证 期 模拟 结果 


Fig. 5 Simulation results during model calibration period and validation period 


未 出 现 较 大 降水 时 ,也 会 出 现 春季 小 洪峰 ,此 时 主 
要 是 由 于 流域 上 游 祁 连 山 东 段 的 山 间 积 雪 融化 导 
致 。 如 2011 年 4 月 20 日 ,此 前 4d 及 当天 未 出 现 降 
水 ,出 现 了 13.8 m?+s 1! 的 流量 峰值 。 由 于 当时 气温 
逐渐 升 高 ,山上 的 积 雪 逐渐 融化 ,形成 了 此 次 春季 
非 降雨 成 因 的 融雪 径流 。 改 进 的 新 安江 模型 在 夏 
季 洪 峰 的 模拟 上 效果 很 好 ,这 可 能 是 因为 在 西 营 河 
流域 , 产 流 层 较 薄 , 多 数 情况 下 当日 降水 会 产生 当 
日 出 流 ,流域 汇流 延 时 较 短 。 在 验证 期 ,总 体 模拟 
效果 普遍 不 如 率 定 期 ,但 春季 径流 仍然 能 够 模拟 较 
好 ,说 明 模型 加 入 的 融雪 模块 运用 较 好 ,以 气温 为 
基础 足以 计算 出 研究 区 的 春季 径流 。2017 年 的 夏 
季 洪 水 模拟 较 好 ,而 2018 年 7 月 .8 月 的 径流 ,虽然 
趋势 一 致 ,但 多 出 现 模拟 值 较 高 的 情况 。 对 照 降 水 
数据 可 以 发 现 , 在 2018 年 修正 的 降水 仅 为 9 月 降 
水 。 这 是 由 于 在 降水 修正 时 ,没有 将 非 较 大 值 的 降 
水 考虑 在 修正 范围 内 ,修正 的 两 个 约束 并 不 能 识别 
较 小 降水 的 趋势 性 问题 ,这 也 是 在 后 续 降 水 数据 使 
用 中 ,需要 深入 探讨 的 问题 。 


3 结论 


<i 


本 文 考虑 到 研究 区 实际 情况 将 新 安江 模型 加 
入 融雪 模块 ,考虑 融雪 对 径流 的 影响 。 男 一 方面 ， 
由 于 研究 区 水 文 站 点 少 ,采用 格 点 降水 和 气温 数据 
集 作 为 改进 的 新 安江 模型 的 驱动 数据 。 运 用 过 程 
中 发 现 异 常 点 并 进行 了 识别 方法 及 修正 方法 的 探 
索 ,将 格 点 降水 进行 修正 作为 模型 驱动 ,进行 水 文 
模拟 ,进而 分 析 新 安江 模型 改进 前 后 、 格 点 降水 数 
据 修 正 前 后 的 模拟 效果 ,得 到 如 下 结论 : 

(1) 在 站 点 资料 缺乏 的 干旱 半 干 旱 区 ,利用 全 
国 范围 内 制作 的 气象 数据 集 解决 由 于 资料 缺乏 , 数 
据 集 存 在 一 定 误差 的 问题 ,运用 降水 -径流 关系 曲 
线 修正 少数 偏差 较 大 的 格 点 降水 后 ,能够 提高 数据 
可 用 性 ,模拟 的 确定 性 系数 在 8 a 年 中 有 6a 均 有 提 
高 ,其余 2a 差 异 不 大 。 

(2) 将 融雪 模块 耦合 新 安江 模型 ,使 用 的 模型 
驱动 数据 较 少 ,运用 简单 。 新 安江 模型 在 西 营 河流 
域 模拟 流量 过 程 的 年 确定 性 系数 均 在 0.5 以 上 ,经 
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Runoff simulation in an arid area using the Xin’ anjiang 


model coupled with snowmelt 
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Abstract: Using datasets of the daily grid-based precipitation and temperature from 2011 to 2018 in the upstream 
research area of Jiutiaoling hydrological station in the Xiying River Basin, we simulated the daily runoff process 
in the study basin. We propose a new method for the identification and correction of abnormal grid precipitation, 
which is based on the relationship between precipitation and runoff. Moreover, we created a combination of the 
snowmelt module and the three water source model of the Xin’ anjiang model, then analyzed the applicability of 
the Xin’ anjiang model coupled with the snowmelt module in the Xiying River Basin. Two evaluation criteria, 
deterministic coefficient and relative error, were used to evaluate the simulation results. The simulation results, 
after precipitation correction and considering snowmelt, were compared with those without precipitation 
correction or snowmelt consideration. The correction of grid precipitation data, based on the relationship between 
precipitation and runoff, improved the accuracy of the runoff simulation, and produced an improved simulation 
certainty coefficient in 75% of the years. The Xin’ anjiang model with a snowmelt module produced a better 
simulation of the Xiying River Basin than that without a snowmelt module in 75% of the years. Additionally, the 
simulated certainty coefficient was above 0.6 for more than 87% of the years. This study can act as a reference 
for the application of the model in the Xiying River Basin. 


Keywords: Xiying River Basin; grid precipitation data; snowmelt module; improved Xin’ anjiang model 


